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We consider macroscopic correlations in a bipartite system consisting of 2N particles described 
by a generalised probabilistic theory. In particular, we discuss a case of N PR-boxes shared between 
two parties. We characterise macroscopic measurements as collective measurements of the same 
property on all the boxes in the same region. Such measurements are assumed to reveal only 
the average value of the measured collective property. We show that for two measurements per 
observer and N >2 there always exist a joint probability distribution explaining all observable data 
and therefore the system admits local hidden variables. Next, we generalise this result to include 
measurement of fluctuations and additional measurement settings. Finally, we discuss our result in 
the context of previous works in which contradictory results were presented. 

PACS numbers: 


I. INTRODUCTION 

Macroscopic systems consist of many particles. It is ex¬ 
tremely difficult, if possible at all, to register individual 
properties of each microscopic constituent. Therefore our 
observations are limited only to collective properties of 
multi-particle ensembles and it is an interesting question 
whether microscopic effects, like strong non-classical cor¬ 
relations, can be still observed at the macroscopic level. 

The above problem has been considered in literature 
before, see for example [IHI2- In particular, it was shown 
that in the macroscopic limit, where measurements are 
limited to reveal only average values [2] and its fluc¬ 
tuations [T], collective quantum correlations become lo¬ 
cal, although individual pairs of particles can manifest 
nonlocality. It is still an open question whether similar 
macroscopic locality [T] also holds for stronger than quan¬ 
tum correlations m- This problem has been recently 
addressed by Rohrlich m and Gisin m who argued 
that such correlations are unphysical since they allow for 
signalling in the macroscopic limit. We show the con¬ 
trary: there exist a local and realistic macroscopic limit 
of stronger than quantum correlations without signalling. 

We generalise our previous result [2], which applies 
only to quantum systems, and show that average val¬ 
ues of macroscopic measurements on a bipartite system 
of 2N particles, whose correlations are described by gen¬ 
eralised probabilistic theories, admit a local and realistic 
description. In particular, we show that for two mea¬ 
surements per observer for A" > 2 there exists a joint 
probability distribution (JPD) for all observables [14] re¬ 
producing measurable average values. Next, we make an 
extension to more measurement settings and to macro¬ 
scopic measurements that reveal information about sta- 
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tistical fluctuations. Finally, we compare our results to 
that of Rorlich and Gisin iniiig. 

II. MACROSCOPIC MEASUREMENTS 

Let us define macroscopic measurement. Following 
previous works HHH], it is a collective measurement of 
the same property on all N particles, i.e., 

A = aW+a(2) + ...+aW, (1) 

where is a microscopic observable on the i-th particle. 
Here we assume that microscopic observables are binary 
with ±1 outcomes. There are 2^ microscopic states of 
such a system and they give A + 1 different outcomes of a 
macroscopic measurement: A, A — 2 ,...,— A. Although 
A +1 ^ 2^ we assume that it is still impossible to trace 
all the macroscopic outcomes in realistic measurements 
and that we are limited to measurements of the A’s aver¬ 
age value (A). Measurements of A’s fluctuations will be 
discussed later. 

To justify our definition of macroscopic measurement 
we use an analogy to macroscopic measurement of a light 
beam’s polarisation. Although the beam is composed of a 
large number of photons, the macroscopic measurement 
of its polarisation (say horizontal and vertical compo¬ 
nents) reveals only intensities Ih and Iv- In such mea¬ 
surements one can also observe fluctuations around the 
mean value but these fluctuations are mainly caused by 
the internal noise of measurement devices rather than by 
actual fluctuations of the observed property. 

An interesting problem turns up when a macroscopic 
system is divided into two parts. One can measure dif¬ 
ferent property on each part, say A and B (defined in 
a similar way), and ask if the correlations between these 
properties admit local hidden variables. It is a non-trivial 
question since we know that the correlations between the 
microscopic constituents in these divisions may not admit 
such a description. 
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There is one important obstacle in an experimental 
implementations of this scenario - there is no known 
method to measure {AB) without invoking nonlocal mea¬ 
surements. In order to measure (AB) one uses a probe 
that interacts with both partitions, say first with A and 
later with B. The final measurement of the probe re¬ 
veals (AB) (a possible implementation technique would 
be similar to the one considered in m)- In this case 
there is an obvious locality loophole, since the probe can 
transfer some information from A to B. On the other 
hand, if one used two separate probes, a different probe 
for each partition, the first probe would reveal (^4) and 
the second (B) but not {AB). 

Because of the locality loophole the above scenario 
would rather test the contextuality [16] of the system 
than its lack of a local and realistic description. However, 
even if local measurements of (AB) were possible, we 
show that there exist a JPD reproducing average values 
of macroscopic measurements. Thus, even if microscopic 
correlations are not local and realistic, the nonclassical 
effects resulting from these correlations are unobservable 
on a macroscopic level - they can be reproduced by a 
system admitting local and realistic/non-contextual de¬ 
scription. 


III. MAIN RESULT 


A. Macroscopic correlations 


Let us consider 2N particles shared by Alice and Bob. 
In general, the correlations between these particles may 
not be classical and not even quantum, but can be de¬ 
scribed by general probabilistic theories m- Each par- 

(k) 

tide can be measured in two different settings Uq ^ and 
for the Alice’s k-th particle and 6 q ^ and for the 
Bob’s 1-th particle. We assume that all measurements 
are binary with ±1 outcomes. 

We can calculate the following probabilities 


_ ^( 1 ) ^( 1 ) - ^( 1 ) - yW) 


( 2 ) 


where ik^ji =0,1 label measurement setting {k^l = 
1,...,A) and = ±1 are corresponding out¬ 

comes. We will write ® as 




pW (1) 

Sat ’ yji 




(3) 


We also assume that these probabilities are non- 
signalling, i.e., the same marginal probability obtained 
from two different global probabilities are equal. For ex¬ 
ample 


However, as macroscopic measurements 0 do not ad¬ 
dress each particle individually, we are limited to a sub¬ 
set of probability distributions in which all settings in 
a partition are the same, i.e., ii = ... = = i and 

jf'i = ... = jAT = j. Thus, we are interested in the follow¬ 
ing probabilities 


p{x. 


( 1 ) 





(5) 


In our scenario there are 
ments 

four macroscopic measure- 

N 

N 

II 

M 

(6) 

k=l 

1 = 1 


where i, j = 0,1. The measurement of {AiBj) yields 


N 

(AB^) = ( 7 ) 

k,l=l 

where 

{afhf) = p{xf'’ = +l-,yf = +1) 

+ p{xf^ = = - 1 ) 

- p{xf^ = +l]yf = -1) 

- p{xf'^ = -l-,yf = +1), (8) 

and p{x^^^ \y^p) is a marginal of (j^. 

An interesting feature of macroscopic measurements 
0 is that the macroscopic outcome stays the same un¬ 
der the particle exchange. In particular, if we consider 
all possible permutations of particles in each region then 
microscopic correlations between every pair will be given 
by an effective correlation 


1 ^ 


k,l=l 

resulting from an effective probability 


N 


(9) 


pixf, yj)eff = = a;*; yf = %)• (lo) 


kd=l 


Thus 


{AiBj) — N {aibj)eff- (11) 

It is clear that {AiBj) is simulated by effective correla¬ 
tions between a single pair of particles, see Fig. 


B. Joint probability distribution for correlations 


E 

' • ’ ’ yji ’ • 

II 


E 

' • ’ ’ yji ’ • 

• • ’ yjN > ■ 

(4) 

x[^^=±l 





Although individual particle addressing is not allowed 
in the macroscopic limit, marginal probabilities of ^ 
exist 




( 12 ) 
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FIG. 1: Effective description of macroscopic measurement of 
{AB). a) The original distribution describing arbitrary cor¬ 
relations. b) Macroscopic measurements cannot distinguish 
between the original distribution and its symmetrised ver¬ 
sion. c) Symmetric correlations are the same for every pair, 
d) Effective simulation with a single pair. 


We can use them to construct a symmetric probability 
distribution 


P{xq^ Xi] yi^gym — 


1 


(13) 


tsym 7V2(A'’-1)2 

kpk^l 

m,n|m7^n 

This distribution naturally arises under the symmetrisa- 
tion over all particles in each region. 

However, it is straightforward to show that 

p{xo;yo)eff = Y P^^o,xi;yo,yi)svm, 

Xi,yi 

p{xo;yi)eff = E] ^ 0 , 2 / 1 ) sym ■) 

Xi,yo 

p{xi;yo)eff = Y P{xo,xi;yo,yi)sym, 

xo,yi 

p{xi;yi)eff = E] P(^o,a;i;yo,2/i)s 


(14) 


xo,yo 


Therefore, (13) is a joint probability distribution that 
reconstructs measurable average values and hence these 
values admit a local and realistic/non-contextual descrip¬ 
tion [14]. 


C. Macroscopic fluctuations 

We now deal with macroscopic fluctuations around the 
average values as in principle they can be measured too. 


The probability distribution (13) recovers macroscopic 
correlations {AiBj) and local averages {Ai) and (Bj) but 
not their fluctuations 


A(A,) = ^(Af) - (15) 

and 

A{AiBj) = (16) 

Therefore, we need to look for a way to reconstruct {A‘f) 
and HABj)^). 

First, note that 


N 


(A?) = ^ {a^aY =N+ Y K ) (17) 

kpkj^l 


k,l=l 


and 


N 


{{AiBjf) = Y 

k,l,m,n=l 

+ iV ^ (afaf^iV ^ {bphf) 

k,l\kj^l m,n\mj^n 

+ Yi ( 18 ) 

kpk^l 

m,n\m^n 


where the 4-point correlations are given by the marginals 
of (j^, i.e., p{x\^\ ). The bipartite correla¬ 
tions and are given by p{x\^\xf^) 

and p{yj ^^, yj ^^) respectively. 

As before, the macroscopic measurements cannot dis¬ 
tinguish between the original probability distribution and 
the effective one obtained via symmetrisation 

p{xi,x'i;yj,y'j)eff = _ 1)2 ^ (1^) 

Y pAP 

kpk^l 

m,n\mj^n 


Primes are used to denote two different outcomes 
on two local observables. Although we distin¬ 
guish between these outcomes, the effective probabil- 
ities obey p{xi,x'i;yj,y'j)eff = p{x'i,Xi-,yj,y'Y)eff = 
p{xi , x '^; 2 /' ,yj)eff= P{x'i , Xi ; y', yj ) e// • 

As a result {{AiBjY) is effectively simulated by 
two pairs of particles admitting the distribution 
p{xi^x'pyj^y'j)eff^ see Fig. [21 This distribution yields 
{aia'^b jb'^) ef f^ {aia'^)eff^ {bj^eff and allows us to ex¬ 
press (17) and (jlSj) as 


(A 2 )=iV(l + ( 7 V-l)(a,a')e//) 


( 20 ) 


and 


{{A,BjY)=N\N-1)x (21) 
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FIG. 2: Scheme representing an effective description of a 
macroscopic measurement of {{ABY). a) The original dis¬ 
tribution describing arbitrary correlations, b) Macroscopic 
measurements cannot distinguish between the original distri¬ 
bution and its symmetrised version, c) Symmetric correla¬ 
tions are the same for every two pairs, d) Effective simulation 
with two pairs. 


Finally, we write a symmetric probability distribution 

((iV-4)!)^ 


p(a;o, a;o, a;i, x'l; yo, 2/0, yi, 2/i) 


sym 


{my 


where denotes the sum over k Y ^ ^^d m Y 

n 7 ^ m' 7 ^ n'. The probabilities p{xi^x[]yj^y'Yeff are 
marginals of (22). Moreover, the probabilities ( [1Q| ) are 
also marginals of (22). This confirms that ( [2^ is a JPD 
that reproduces average values and their fluctuations. 

At this point we need to comment on an important 
property of our approach that differs from the standard 
one m- For average values we showed that the origi¬ 
nal system can be effectively simulated by a single pair 
of particles whose correlations are described by the JPD 
(13). To obtain the original average value from the simu¬ 


lated one we only need to rescale the simulated outcome 
by the factor of N^. Therefore, the average values stem 
directly from the JPD. 

On the other hand, to simulate fluctuations one needs 
two effective pairs described by the JPD (22). However, 


the original fluctuations cannot be obtained by a simple 
rescaling of fluctuations of the simulating system. The 
original fluctuations do not stem directly from the JPD, 
but are its functions (^ and ( [^ . From the practical 
point of view the local and realistic/non-contextual sys¬ 
tem of four particles, together with two local devices com¬ 
puting functions of the JPD, can be effectively used to 
simulate averages and fluctuations of macroscopic mea¬ 
surements on a bipartite system consisting of many non- 
classically correlated particles. 


D. Higher moments and more measurement setups 

The above results can be generalised to explain the be¬ 
haviour of higher moments {{AiBjY) and scenarios with 
more than two measurement setups per observer. Also, 
note that the above discussion does not depend on the 
number of measurement outcomes in a given setup. 

Without going into details, let us note that our method 
relies on the symmetrisation of the probability distri¬ 
bution and on choosing a proper subset of particles 
on which we measure all the observables. For exam¬ 
ple, for /c = 3 one needs to measure correlations of 
the form This means that in 

case of two measurement setups per observer one needs 
3 + 3 + 3 + 3 = 12 particles to obtain the corresponding 
JPD. 

In general, in order to obtain a JPD simulating the 
k-th moment in case of sa and sb measurement setups 
one needs k{sA sb) particles. This shows that in or¬ 
der to look for macroscopic nonclassicality one needs to 
make very accurate measurements of fluctuations to re¬ 
veal higher moments k{sA + 5^) > N. Indeed, if one 
could measure sufficiently high moments, one would gain 
enough insight into a macroscopic probability distribu¬ 
tion p{Ai = Xi]Bj = Yj)^ where Xi^Yj = N^N — 
2,..., — A", and would be able to apply methods of Ref. 
[3] to observe violation of local realism. This seems to be 
unrealistic in a macroscopic limit with Avogadro number 
N ^ 10^^ of particles. 


IV. EXAMPLE - N PAIRS OF PR-BOXES 

Here we show how to apply our result to a specific 
scenario. Let us consider a bipartite system consisting of 
N pairs of PR-boxes. PR-box (Popescu-Rohrlich-box) is 
a hypothetic bipartite system manifesting stronger than 
quantum correlations without signalling m- 


A. Correlations 

Each PR-box pair is described by correlations 

{aibj)pR = (23) 

and the corresponding probabilities 

+ ( 24 ) 

Local marginal probabilities are maximally random 

PPR{xi) = PPRiVj) = hence (oj) = (bj) = 0. In the 
case of N pairs each PR-box pair is independent, i.e., 
probabilities factorise. For example, for any two pairs k 
and I we have 

+ = PPR{Xi\ym)PPR{Xj;yn)- (25) 
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The above implies that there are no correlations between 
boxes corresponding to different pairs since 

= PPR{Xi)ppR{yn) = (26) 

therefore 

(af = 0. (27) 

Correlations between macroscopic observables are 
given by 


N 

{AiBj)= Y,{af^hf'^)=N{aibj)pR = N{-iy-^. (28) 


k,l=l 


We also have {A^) = {Bj) = N. Again, this is because 
= 0 for /c 7 ^ / and m ^ n, whereas 


{af^afy = {A A) 


= 1 . 


The above expressions can be also recovered from a 
symmetric system of four boxes described by the effec¬ 
tive probability distribution Pe//(^i 5 ^^5 %5 ^ 7)5 see (19). 
For all settings these probabilities have only four possTble 
values 


a 


7 


N{N - 1) + 2 (AT - 2){N - 3) 

16Ar(A^-l) ’ ^ ~ 16A^(A^ - 1) ’ 

N{N + ?,)-2 ^ (jV + i)(Ar-2) 

16A^(A^-1) ’ 16N{N-1) ' ^ ’ 


However, we know from the previous section that this For the settings AqBq^ AqBi and AiBq these values cor- 
can be explained by a single effective pair of boxes respond to the following outcomes: 


{AiBj) = N^{aibj)eff = V 


2 / (-1)"' 


N 


(29) 


with the corresponding effective probability distribution 
1 ^ 

Peff{Xi-,yj) = ^ 


k,l=l 




4A^2 


fe=l 

N - l+4ppp{xi;yj) 


1 

- 4 + 


m 

\xi + {-iy-^yj\-l 

4N 


(30) 


B. Fluctuations 


Next, let us derive {{AiBjf). From (18) we have 


{{AiB^f)=N^+ Y, (31) 


k,l\k^l 

m,n\m^n 


This is because a^P) = {bp^^bp^) = 0 for /c 7 ^ / and 

m ^ n. Morevoer, because and b^ are not correlated 
for /c 7 ^ / one has 


{(ABjf) = N^+ Y 

k,l\kj^l 

+ Y ^>7)• 

k,l\k^l 


(32) 


Factorisation of distinct PR-box pairs 
gives {ap\Pbp\P) = {ap^aPb^bf^ = 
{ap^bpp{aPbP) = {ciibj)‘^j^ and therefore 

\2\ _ 7\r2 


{{AB.Y) = N^^2N{N-l){aAjypj, 

= 3N‘^ - 2N. (33) 


• (+5 +: ~): (+5 +): (~: +; +? “)? (+: +); 

P • (+5 +; ~): (~5 +: +); 

7 • (+5 +; +: +): (~5 ~); 

5 : remaining outcomes, 

whereas for the setting AiBi one has 


^ • (+5 +5 “): (+5 +)? (“: +; +5 —)^ (+? +); 

P : (+, +; +, +), (—, —; —, —); 

7 • (+5 +; ~): (~5 +: +); 

5 : remaining outcomes. 


It is straightforward to verify that the above probabilities 
recover (30) as marginals. 


We note that these effective probabilities lead to 
{aia[bjb'j)eff = 1 - 166 = 




(35) 


If we rewrite (pTl) as 


{{AiB^f) =N^ + N\N - l)7a,a'&,-6; )e//, (36) 

we obtain the same result as in 

The effective distribution also recovers local ffuc- 
tuations (A?) given by (20). This is because 
marginals Peff{xi,x[) lead to' 7 ^^a')e// = 0 (similar for 
PeffiVj^Vj))' This proves that our system is effectively 
described by four symmetrised boxes. 


C. Joint probability distribution 


There exist JPDs (13) and ( [^ for which (30) and (34) 
are marginal distributions. Using simple, yet tedious, 
combinatorial techniques we calculated the JPD (22) for 
a system of N PR-boxes and confirmed that it recovers 
measurable averages and fluctuations. It consists of 256 
probabilities that are represented by the ratio of third 
order polynomials in N, therefore for obvious reasons we 





















6 


do not present them here. However, we present the JPD 
for averages (13) that consists only of 16 probabilities. 

The 16 probabilities of the JPD p(xo, xi; ^i) fake 
only two possible values 




N±2 

16 A^ ■ 


(37) 


The value co’+ corresponds to the following events 


(+: +; +5 +)? (+5 +; +: (+5 +: +)? (+? +)? 


Without going into details, let us recall the gist of 
Rohrlich’s reasoning. For a single PR-box if Alice mea¬ 
sures ao then if Bob measured they would observe 
yo = xo, but if he measured hi they would observe 
yi = xq. If it were possible to measure both 6o and 
hi one would have yn = yi. However, if Alice measures 
ai then yo = xi, but yi = —xi, and the possibility of 
joint measurement would show that yo = —yi. 

This microscopic feature seems to affect the variance 
of the macroscopic measurement 

Bo + B, = + 6^) + ... + ^ (40) 


and the value uj- correspond to the remaining ones. From 
this we see that for all settings, except i = j = 1 the 
marginals p{xi;yj) are 


The average value {Bq + Hi) = 0 is independent of 
whether Alice measured Aq or Ai, however this does not 
happen for the variance since 


p(+; +) — p( —; —) — — 

p(+; —) = p( —; +) = 3 lj - + = 

For the setting i = j = 1 one gets 

+) = ~) = + <^+ = 

p(+; —) = p( —; +) = 3cj+ + = 


A^ + 1 

AN ’ 

1 

~Jn~‘ 


N-l 

AN ’ 
A^ + 1 

4A^ ' 


(38) 


(39) 


This reproduces (30) and therefore confirms our claim. 
Note, that for A" > 2 the JPD is a valid probability 
distribution which for A = 1 becomes negative. 


V. DISCUSSION 

Finally, let us discuss our results in the context of 
the previous works. In particular, we focus on the re¬ 
cent works by Rohrlich m and Gisin m who argued 
that stronger than quantum correlations do not admit a 
macroscopic limit. 


A. Previous results 

The argument of Rohrlich is based on an observation 
that although correlations in a single PR-box pair do not 
allow for signalling, a collection of them would signal in 
a macroscopic limit in which a joint measurement of in¬ 
compatible observables is possible. He considers macro¬ 
scopic measurements on A pairs of PR-boxes - the situ¬ 
ation discussed in this paper. Next, he assumes that in 
the classical limit macroscopic measurements allow for 
a joint measurement of observables that are microscop¬ 
ically incompatible. In particular Bq and Bi should be 
jointly measurable. However, his treatment of macro¬ 
scopic measurements shows that, because of strong non- 
classical correlations, the joint measurement of Bq and 
Bi allows for signalling from Alice to Bob. 


{iBo + Bjf)Ao=0{N), ((Bo + Bi)')^, =0. (41) 

Therefore, the signalling can occur via fluctuations. Be¬ 
cause of that Rohrlich claims that macroscopic measure¬ 
ments on PR-boxes do not admit a classical limit and 
hence are unphysical. 

Interestingly, Rohrlich also shows that even quantum 
correlations lead to a signalling unless one assumes that 
Bq and Bi are not scalars. Such apparent signalling in 
quantum theory is in fact analogous to the (in)famous 
von Neumann ” silly mistake” m on which he founded 
his proof of the impossibility of hidden variables - the 
outcome of a measurement of Hq + Hi cannot be assumed 
to be a sum of the corresponding outcomes. What hap¬ 
pens is that in quantum theory expressions like Bq-\-Bi or 
BqBi are completely new objects. However, this solution 
of the quantum signalling paradox weakens Rohrlich’s 
main claim, since there is no reason why for PR-boxes 
a sum or a product of Bq and Bi should not be a new 
object either. 

Another argument was used by Gisin [12] to show that 
there is no macroscopic limit of stronger than quantum 
correlations. Gisin first uses quantum formalism to study 
weak measurements on A singlet pairs. Weak measure¬ 
ments allow him for a joint measurement of two incom¬ 
patible observables. He shows that a signalling in quan¬ 
tum case is prevented by the unavoidable fundamental 
noise in the measurement process. Next, he considers 
noisy isotropic PR-boxes and simply assumes that in a 
macroscopic limit a JPD pg{Aq, Ai, Bq^ Bi) exists and 
that it is Gaussian because of the central limit theorem. 
Finally, using the argument that Gaussian distribution is 
nonnegative if and only if the corresponding correlation 
matrix is nonnegative 

/{AqAq) {AqAi) (AqBq) (AoHi)\ 

{AiAq) {AiAi) {AiBq) {AiBi) 1 . . 

(BqAq) (BqAi) (BqBq) (BqBi) ’ 

\{BiAq) {BiAi) {BiBq) {BiBi)J 

he shows that the above matrix is nonnegative when the 
correlations between the PR-boxes are at most as strong 
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as quantum ones. Gisin managed to obtain this result 
without estimating (AqAi) and (^o^i) because these 
values are bounded by the other measurable values. Let 
us also remark that for any probability distribution the 
covariance matrix has to be nonnegative [T], but in the 
case of PR-boxes local averages are zero and the covari¬ 
ance matrix is equivalent to the correlation matrix. 

B. Why our result contradicts Rohrlich’s and 
Gisin’s ones 

Let us first summarise how we arrived at our result. 
Our approach consists of the following major steps: 

• Assumption: macroscopically observable data is 
limited to average values (later this assumption is 
relaxed to include variances and higher moments 

• Observation 1: from the point of view of macro¬ 
scopic measurements the same results would be ob¬ 
tained if one symmetrized over all particles in each 
region. 

• Observation 2: for symmetric systems measurable 
data can be described by a small subset of parti¬ 
cles for which one can find a JPD. This small local 
and realistic/non-contextual subset can effectively 
simulate macroscopic measurements on the whole 
system. 


allow to recover them. It originates from an averaging 
over permutations of all particles in the system, hence 
its derivation is an irreversible process in which some in¬ 
formation is lost. Moreover, we showed that our JPD 
approach brakes down if too much information about 
microscopic properties is available, since it only allows 
to recover limited number of moments {{AiBj)^). Never¬ 
theless, because of this loss of information new properties 
like Bq + Bi can emerge and the original system becomes 
simulable by a much smaller local and realistic/non- 
contextual system. In addition, the JPD is compatible 
with Bq and Bi derived from microscopic properties, but 
it is not compatible with Bq + Bi derived by Rohrilch, 
hence our extension of PR-box formalism is different than 
his. 

Next, we compare our result with Gisin’s one. If one 
evaluates the correlation matrix using our methods one 
gets the diagonal elements {AiAi) = (BjBj) = N and 
the off-diagonal elements {AiBj) = A^(—These val¬ 
ues were directly calculated from the microscopic prop¬ 
erties and were shown to be recoverable from our JPD. 
The remaining off-diagonal elements {A^Ai) and (^o^i) 
cannot be calculated from the microscopic properties and 
therefore they need to stem from the JPD. If one defines 
these properties similarly to {AiBj)^ i.e., 

(AqAi) = N‘^ {aoai) ef f (43) 

and 


The above leads to a macroscopic limit of non-classical 
correlations that does not allow for signalling and that 
includes stronger than quantum correlations. 

We begin with the divergence of our and Rohrlich’s 
results. Because there is no reason why observables in 
generalised probabilistic theories should not be repre¬ 
sented by non-scalar entities and since the PR-box for¬ 
malism defines only what happens if local microscopic 
measurements are either bo or bi , one has a lot of freedom 
while making assumptions about the nature of a joint 
measurement of two macroscopic properties Bq and Bi. 
Therefore, there could be many extensions of the under¬ 
lying formalism that include new elements like Bq Bi. 
Each extension can be based on different assumptions. 
Rohrlich’s claim is based on an assumption that for PR- 
boxes the values of Bq + Bi and BqBi are determined by 
values of microscopic properties and b ^^. His partic¬ 
ular extension is later shown to be unphysical, however 
this does not mean that every extension must be unphys¬ 
ical. 

On the other hand, in our approach and b^^ only 
define Bq and Hi, but not Bq + Bi and BqBi. The 
latter macroscopic properties can be derived from the 
JPD which in our case is more fundamental. What is 
important, although the JPD is calculated form the mi¬ 
croscopic properties of the original system, it does not 


{BoBi)=N^{boh)eff, (44) 


where the effective two-particle correlations are derived 
from the marginals p{xo,xi) and p{yo^yi) of the JPD 
(22), one finds that {AqAi) = {BoBi) = 0. However, 
having all entries of the correlation matrix one finds that 
it has two negative eigenvalues ^"(1 — v^), therefore we 
arrive at a paradox, since this matrix cannot correspond 
to a valid probability distribution. 


This paradox is resolved when we realise that in our 
case the actual entries of the correlation matrix do not 
stem directly from the JPD, but rather from its func¬ 
tions, and for different entries we use different functions. 
If these entries were calculated exactly from the JPD 
the matrix would be nonnegative. Once again we stress 
that in our approach measurable data is not directly re¬ 
coverable from the JPD, as is assumed by Gisin. The 
main point is that we show that there exist a local and 
realistic/non-contextual system (described by the JPD) 
that is capable of simulating the original system. 
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